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We consider the pion structure in the region of low and moderately high momentum transfers: at low Q 2 , 
the pion is treated as a composite system of constituent quarks; at moderately high momentum transfers, 
Q 2 = 10-^25 GeV 2 , the pion form factor is calculated within perturbative QCD taking into account one-gluon 
hard exchange. Using the data on pion form factor at Q 2 < 3 GeV 2 and pion axial-vector decay constant, 
we reconstruct the pion wave function in the soft and intermediate regions. This very wave function combined 
with one-gluon hard scattering amplitude allows a calculation of the pion form factor in the hard region 
' Q 2 = 10 -T- 25 GeV 2 . A specific feature of the reconstructed pion wave function is a quasi-zone character of the 
gg-excitations. On the basis of the obtained pion wave function and the data on deep inelastic scattering off 
the pion, the valence quark distribution in a constituent quark is determined. 

o 

1 Introduction 

Perturbative QCD gives rigorous predictions for exclusive amplitudes, in particular form factors at asymptoti- 
£f) \ cally large values of Q 2 Q . For the pion form factor defined as 

<7r(P')\J p \7r(P)>=(P' + P) p FjQ 2 ) (1) 

co ; 

, the pQCD result takes the form 

in 

FlT (Q 2 ) = 8nf 2 ^l ll+ V C N (ln9l) l+OCI/O 1 ). 12) 
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O || Here 77V are the known positive numbers calculated within pQCD, fj, is a constant about 1 GeV dividing the 
perturbative and nonperturbative regions; f n — 130 MeV is the pion axial-vector decay constant, and the 
Cn are expressed through the soft-region nonperturbative wave function of the pion. In the series (||), the 
contribution of diagrams with internal lines having virtualities above /z are taken into account perturbatively, 
while all the exchanges with lower virtualities are absorbed into the set of soft nonperturbative wave functions of 
the pion Fock components. The terms of the order 1/Q 2 come only from the valence quark-antiquark component 
of the pion Fock state, whereas other components give the terms 0(1/Q 4 ). The series (Q) involves both the 
leading and subleading logarithms. In the leading logaithmic approximation (LLA) the expression (|^) can be 
rewritten in the form 



1 

F„(Q 2 )= I dxdx'^(x,Q 2 )T B (x,x',Q 2 )^(x',Q 2 ) (3) 





where (j>(x,Q 2 ) is the leading twist wave function (distribution amplitude) which describes the longitudinal 
momentum distribution of valence quark-antiquark pair whose relative transverse momentum is less than Q, 
and 

W ,g.>^*a£gp (4) 

is the amplitude of the hard interaction of the two free quarks in the Born approximation (Fig.^a) . 

The distribution amplitude at large Q is related to the soft pion distribution amplitude (f) v (x, /j, 2 ) by the 
gluon ladder evolution kernel K(x,y,Q 2 / fi 2 ) as follows 

fa{x, Q 2 ) - / K(x, y, Q 2 /fj, 2 )My, ^) dy (5) 
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Figure 1: a. Hard scattering amplitude for two free quarks in the Born approximation, b. Hard scattering 
process in the LLA. c. Soft form factor. 



The soft distribution amplitude is connected with the pion axial vector decay constant via the relation 

j fafa M 2 ) das = U, f„= 130 MeV (6) 

So, in the LLA hard scattering off the pion described by the expressions (||) and (||) (see Fig.|l]b) has a clear 
physical interpretation within the hard scattering picture ||, namely: The initial pion transforms into quark- 
antiquark pair with a small relative transverse momentum < /i and the longitudinal momentum fractions x and 
1 — x, respectively: this stage is described by 4>(y,fj, 2 ). Next, the quarks are coming closer to each other to 
the distances 1/Q and increase their virtualities via ladder gluon exchanges, described by K. Then, the hard 
interaction of almost free quarks with the external current occurs (Tg), and inverse evolution to low virtualities 
(K) with a subsequent pion formation {<f>(y' , H 2 ))- A rigorous QCD result is that any soft distribution amplitude 
evolves at large Q 2 to the universal function 

C s (z) = M*, Q 2 °°) = 6 f* x (! - ( 7 ) 

Substituted into @, this function gives the leading term in (|^). 

Unfortunately, this beautiful picture does not work at momentum transfers accessiblle to present-day ex- 
periments: the momentum transfers are not large enough. Trying to adjust this perturbative calculations for 
moderately high momentum transfers, Chernyak and Zhitnitsky || assumed the logarithmic and power correc- 
tions (including a purely soft contribution shown if Fig.|l|) to the Born term to be not essential, and the pion 
form factor to be described by 

F^(Q 2 )= f ^^(x,/i 2 )T B (x,x',Q 2 )dT'^(a;',Ai 2 ) (8) 
Jo 

with Tb still given by (|j), but some modified distribution amplitude cj)^. Describing the available pion form 
factor data by the formula (JsJ) , they came to the soft distribution amplitude of the form 

^(x, ii w O.bGeV) = 30/^(1 -x)(l- 2x) 2 (9) 

However, arguments against such an approach were put forward by Isgur and Llewellin Smith B. The 
problem is that the wave function of the form (|^) strongly emphasizes the end-point region of x ~ which 
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Figure 2: The expansion of the pion form factor in the series over g 2 . 

was estimated to give 70 ~ 90% of pion form factor at Q 2 — 10 GeV 2 . Recall that only the exchanges with 
virtualities above fi 2 were considered perturbatively, otherwise the corresponding subprocesses were referred to 
the soft wave function. In the end-point region at moderately high Q 2 , the gluon virtuality xx'Q 2 is not large 
enough to justify the perturbative treatment. Hence, the end-point contribution should be rather referred to 
the nonperturbative one. The large contribution coming from small x means that in fact the wave function 
with strongly emphasized end points has picked up a good portion of the nonperturbative contribution. So the 
latter turns out to be not small that contradicts to the inital assumption. 

The arguments of 0] were supported by by recent applications of QCD sum rules [f|,||: the end-point con- 
tribution remains numerically important at least up to Q 2 « 10 GeV 2 , although parametrically it is suppressed 
by an extra power of 1 /Q 2 . The problem of a correct extraction of the end-point contribution to hard scattering 
amplitudes was studied by Li and Sterman H and Radyushkin Their results give additional arguments 
against the application of the strategy of ref. to hadron form factors at intermediate momentum transfers. 

In addition, some problems are encountered when wave functions with an emphasized end-point region are 
applied to deep inelastic scattering. It was discussed || that the valence quark contribution to the deep inelastic 
structure function calculated with such distribution amplitudes exceeds the available data at large x, whereas 
there should be a room for other Fock components. 

The attempts to describe the pion form factor in the region Q 2 = 3 -j- 10 GeV 2 taking into account only the 
perturbative hard scattering mechanism were not successful. The nonperturbative contribution in this region 
is obviously not small. Our goal is to consider the pion form factor at moderately large Q 2 allowing for both 
the perturbative and nonperturbative contributions starting with low Q 2 and advancing to higher values. 

Investigations of soft hadron processes in past decades have demonstrated constituent quarks to be relevant 
objects for describing hadron structure p0[pl]| . 

The nonperturbative contribution to the pion form factor was considered within the framework of the QCD- 
inspired constituent quark model (|[2]][jl3| |l4j [[l5j |l6| ) , the pion form factor being represented by the diagram 
of Fig|l]c. 

In the approach suggested here, we move to the region of large Q 2 starting with small values. That is, we 
take into account the diagram of Figj^a and, as the following step, the diagrams of Figs.||b,c. The diagrams 
of Fig.pjb and Figj^c correspond to the terms with the minimal N = in the series (|J). In fact, we recast the 
series O) for the pion form factor as a series over a s . Such an expansion is relevant at moderately large Q 2 . 

To be more quantitative, the procedure is as follows. The expansion for the pion form factor in a series over 
a s reads 

J» = [q(0h p q{0)]^ + a *J <fei<fe2[«(0)8(*i)] / ,»|0 > r!(*i,js a ) < 0\[q(z 2 )q(0)]^ + 0(a 2 s ) (10) 

To make this expansion meaningfull, the operator product expansion was performed, i.e. all the field operators 
were decomposed into soft and hard components. The subscript implies that the subdiagrams for corre- 
sponding operators involve only lines with virtualities below fj, 2 . The contribution of the region of virtualities 
below /i 2 is described by the soft wave functions, whereas the contribution of larger virtualities is represented 
by the hard scattering block Tg Q . 

Let us denote the corresponding contributions of the first and second terms in the r.h.s. of ( |l0| ) as F ss 
(Fig.|a) and 2F sh (Figs.|b,c), respectively. Then one finds 

F-n = F ss + 2F sh + 0(a 2 s ). (11) 

The last series actually corresponds to dividing the pion light-cone wave function \I/ into two parts such that 
^ soft is large at s = x 7^t x \ < s o while ^hard prevails at s > sq. We perform such a decomposition of the wave 
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function using the simplest ansatz with the step-function: 

* = ^softO(s - S )+ ^ h ardO(s - S ) (12) 

According to jl0|), ^hard is represented as a convolution of the one-gluon exchange kernel V as with ^ soft 

*kri = ^®f so/i . (13) 

The soft-soft contribution F ss in ([ll]) is the usual quantity calculated within constituent quark models, 
whereas the soft-hard term F sh relates to the one given by the hard scattering mechanism. The soft-soft 
contribution includes the Sudakov form factor of the quark. 
We obtain the following results: 

(i) The pion form factor calculated in the region of Q 2 from to 20 GeV 2 describes well the available data (Fig.||). 
The soft -soft contribution is found to give more than a half of the form factor in the region Q 2 < 20 GeV 2 
and is not negligible up to Q 2 < 30 GeV 2 . However, the particular numbers depend on how we define the 
boundary of the soft and the hard regions. We assume an extended soft region for s < 9 GeV 2 , that yields a 
large contribution of the soft-soft form factor until very high momentum transfers. 

The transverse motion in the soft-hard term turns out to be important. At the same time, the Sudakov 
suppression is not large in the kinematical region of momentum transfers where the soft-soft term dominates. 

(ii) The soft pion wave function which has been a variational quantity of our consideration is found to have a 
quasi-zone structure: it is large at low s < 2 GeV 2 , then it almost vanishes at 2 GeV 2 < s < 4.5 GeV 2 , and 
has a bump at 4.5 GeV 2 <s<9 GeV 2 (Fig.|). 

(iii) The pion structure function is expressed through pion soft wave function and constituent quark structure 
function. By describing the data on pion valence quark x-distribution, we find the parameters of the x- 
distribution of valence quark inside a constituent quark to be in a qualitative agreement with Reggeized QCD- 
gluon intercept calculation . 

The paper is organized as follows: 
The Section 2 considers the pion as quark-antiquark bound state within the light-cone technique jl?]] refor- 
mulated as dispersion relation integrals and presents the expressions for the pion form factor and and quark 
distribution in deep inelastic scattering. All necessary technical details relevant to pion description are given 
in the Appendices. The results are discussed in the Section 3. A brief summary and outlook are given in the 
Conclusion. 



2 Pion form factor and structure function 

The light-cone technique expressed in the form of the deispersion relation integrals |l8[ ] allows constructing 
relativistic and gauge invariant amplitude of the interaction of a composite system with an external vector 
field starting with low-energy constituent scattering amplitude (see the Appendix A). Two-particle s-channcl 
interactions are consistently taken into account both in the constituent scattering amplitude and the amplitude 
of interaction with an external field. In the case of a bound state, its form factor and structure function are 
expressed through form factor and structure function of mass-shell constituents and the vertex G of constituent- 
bound state transition. This vertex is defined by the two-particle irreducible block of the constituent scattering 
amplitude. On the one hand, the dispersion integral representation turns out to be equivalent to the Bethe- 
Salpeter treatment with a separable kernel of a special form, the vertex G being connected with the amputated 
Bethe-Salpeter wave function of the bound state. On the other hand, this approach can be formulated as a 
light-cone description of a bound state with the special form of spin transformation (the Melosh rotation). 
Because of the relativistic invariance, the dispersion integral approach does not face the problem of choosing 
appropriate component of the current for form factor calculation. The function G determines the bound- 
state light-cone wave function pq| . Note also that only amplitudes for on-shcll constituents contribute to 
corresponding amplitudes of the bound state. This guarantees gauge invariance of the derived expressions and 
escapes the problem of constituent amplitudes off the mass shell. All relevant details can be found in the 
Appendix A. 

Our position on the pion structure completely coincides with the viewpoint formulated by Weinberg 
Successes of the quark model allow one to treat quarks as usual massive hadrons with the only difference that 
quarks are subject to color forces which become essential at large distances and keep quarks confined in hadrons; 
in all other aspects these forces are weak and quarks can be treated as real particles. 

In the soft region, quark structure of a pion is described by the vertex 

Q a {P-k)i l5 Q a (k)^ . d2 



G V {P Z ) (14) 
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with a a color index, N c = 3 the number of quark colors, k 2 = to 2 , (P — k) 2 = to 2 , but P 2 = s ^ to 2 . Here to 
is a constituent quark mass. We consider ir + and omit the flavor which gives the unity factor, (the appendix B 
presents a detailed consideration) The soft vertex G v is supposed to be nonzero at s < sq in accordance with 
(|l2|). Once the vertex is fixed, we can proceed with form factor calculation. 

2.1 Soft— soft contribution to pion form factor 

The double dispersion relation integral for a soft-soft contribution to the pion form factor is given by the 
following expression 

Here f c is a constituent form factor, / c (0) = 1. The quantity A w is defined as follows 

— - ( dkidk^dfaSik 2 - m 2 )5(k' 2 ~ m 2 )5{kl - m 2 )5{P - ki - k 2 )S{P' - k[ - k 2 ) 

»7T J 

x Sp ((fcj + m)7 M (Ai + TO)«7s(m - fe^Ts) = 2P^(q)A 7r (s' , s,q 2 ) (16) 



with 

Pu(q) = (P- 

q 

The trace reads 



PM = (P-^r q)„ p 2 = s, p' 2 = s', (p' pf = q 2 . 



^Sp (ik[ + TO)7 M (fci + to)7 5 (to - ^2)75) = m 2 (ki + k' x )^ + /ci M (fc' 1 fc 2 ) + fc^(fcifc 2 ) - k 2ll {qk x ) (17) 

To reveal the relationship between the dispersion integral ( |l5|) and the lihgt-cone technique, we introduce the 
light-cone variables 

k- = -i=(fco - k z ); k+ = -j=(ko + k z ); k 2 = 2k + k- — k\\ (18) 
into the representation (|l^) and use the reference frame 

q+=0, P±=0, q 2 = -q 2 ± <0 : (19) 

The form factor takes the form 

F ss (ql) = - I ' dxcPk x il>(x,kx)il>{x,k 1 . - x qi _)p(x,k ± , qi _)f c (ql), (20) 



7T „ 

where the soft radial light-cone wave function of a pion is introduced 



,, , s G v (s)y/s m 2 + k 2 ± 

w(x,k±) = — — , s = — =r (21) 

V ' irV8(s-ml)^/x{l-x) x{l - x) v 1 

to 2 + k 2 _ — xkxqx 

P 



v/ to 2 + k 2 ^ \J m 2 + (kx — xq±) 2 

The quantity (3 accounts for the contribution of spins. It is different from unity at qx 7^ because both 
the nonspin-flip and spin-flip amplitudes of the interacting quark contribute. The relation F w (0) = 1 is the 
normalization condition for the soft radial wave function 

dx dk\ \il>{x,kx)\ 2 = 1. (22) 

The pion axial-vector decay constant / T related to the tt — > \w decay is given by 



9A 7^J dxdkl 



\J to 2 + k 



i>{x,kx)=U (23) 



where qa is the constituent quark axial- vector coupling constant. From the analysis of the neutron /3-decay 
within various models, the value of qa was found to be in the range from 0.75 (nonrelativistic constituent quark 
model without configuration mixing in the nucleon wave function, to ss 0.33 GeV) up to 1.0 (relativistic quark 
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model with light constituent quarks, m ~ 0.25 GeV |lq]. ) Here we use the relation (|23|) to fix the value of gA 
related to a particular pion wave function such that (23) reproduces the observed value = 130 MeV. The 
value of gA is found to lie in the range gA ~ 0.75 -j- 1 (see the Table |l|). 

The same expressions for the form factor and pion electroweak constant as (p(i|)-(p3"|) were derived in refs 
|L4j-||l(|. In contrast to the mentioned papers where such formulas were applied to describing the total pion 
form factor , we use them for calculating the soft-soft contribution. 

The soft-soft form factor involves the constituent quark form factor which should satisfy the conditions 
/c(0) = 1 and f c (Q 2 ) — > S(Q 2 ) at large Q 2 — —q 2 . Here S(Q 2 ) is the Sudakov form factor which is taken in 
the form || 

Sm . exp (-?f!><W (I)) , O r - ?%=± = J, Q„ - 0.85 + 1 GeV ,24, 
with a s the coupling constant. At low Q 2 we assume a s to be frozen at 1 GeV 2 , namely we set 
«s(Q 2 ) = -^log -1 , A = 0.22 GeV, Q 2 < 1 GeV 2 



9 b V A2 , 

const, Q 2 < 1 Gel/ 2 (25) 



The constituent quark form factor is taken as 



f c (Q 2 ) = 1, Q < Q (26) 
= S(Q 2 ), Q > Q , S(Q 2 ) = 1 



Note that the Sudakov suppression is absent in the soft-hard term. 



2.2 Soft— hard contribution to the pion form factor 

In accordance with ([h]), the soft-hard contribution is described by the two graphs of Fig.|^b,c with one-gluon 
exchanges. The corresponding dispersion relation integral reads 

h _ f ds G v (s) 9( s < So )ds"0( s ">s ) ^ ds'G v (s>) d(s'<s ) 

J ir(s — m%.) 7r(s" — m*) ir(s' — m*) 

In this expression, D(s, s", s' , q 2 ) is a spectral density which takes into account the one-gluon exchange at large 
s" ; the Appendix C presents the details of the soft-hard form factor F sh calculation. Notice that the dispersion 
expression involves only two-particle singularities and neglects three-particle intermediate state related to 
cutting the gluon line. The final result has the form 

™hrJi\ C F f dxd 2 kj_ 6(s" > s ) , . 

dx'd 2 k\ ,, , ,, , , , , , „ 9 a s (-(k-k') 2 ) 
i ? i>{x', k' ± - x'q ± )d(s' < s )Tr{s, s' , s", q 2 )-^ LL 



where 



, i ,. 2 _ (k'± x ~ k±x') 2 + m 2 (x - x') 2 
— (fe — k) — , 

xx 

Tr{s, s', s", q 2 ) = 2s' (a(s" , s, q 2 ){s" + s - q 2 ) + q 2 ) - 4m 2 (a(s" , s, q 2 )(s' + a - q 2 ) + q 2 )) ; 

-<? 2 (s" + s - q 2 ) 



II 2,\ 

a(s ,s,q ) 



(s"-s) 2 -2g 2 (s" + s)+<j 4 



and 



5 „ _ rn 2 + (fcj_ - X9_l) 2 ^ _ m 2 + - .t'^jJ 2 



x(l — x) ' x(l — x) x'(l — x') 

In (|28l ) me is the gluon mass which depends on the gluon momentum squared —{k' — k) 2 and is normalized to 
be of the order of 1 GeV in the soft region | fl9|| . 
One can easily see that F sh (0) = 0. 
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Let us turn to the region of large Q 2 and compare the soft-hard term which dominates the form factor with 
the standard pQCD expression (Q) and (ji|). Using the relations 

Tr -> Am 2 Q 2 x and - (k' 2 - k 2 ) 2 -> xx'Q 2 

we find 

2 1 f 8ir a s (xx'Q 2 ) 



f t(Q 2 )—>— (j>(x,s ) dx— — — —2 — <j)(x', s ) dx' (29) 
t\i j y xx \a£ 

Here we introduced the distribution amplitude <fi(x, sq) normalized by the standard condition 

J 4>(x,s ) dx = f„ (30) 
This distribution amplitude is related to the soft radial wave function (El]) via the relation 



<t>(x,s ) = y-^- / dkl iP(x,k ± )0 — i < so 31 

V27T J y / m 2 + k 2 ± \x(l-X) J 

The expression ( pi] ) differs from (p]) by the factor 1/k 2 , where 

« = <fl(0) = 0,75- 1 (32) 

(see the eq.([78|) of the Appendix B) is the axial-vector decay constant at the level of a constituent quark. If we 
identify constituent quarks with a bare pQCD quarks, then = 1. Here we assume such an identification to 
be not well-justified in the region Q 2 < 10 GeV 2 and hence we allow n to be in the range from 0.75 to 1. 

Actually, a constituent quark is a Fock state involving the components q va i, qvalQQi QvalQQ9t ■ • ■ At asymp- 
totically large momentum transfers only the valence component survives whereas all other components do not 
contribute to the form factor. The problem of the typical Q at which the transfer to the asymptotic regime 
occurs is still opened. It is quite natural that the expression (|2^) gives a larger value for the form factor at 
K < 1 since it includes the contribution of all higher Fock components of the constituent quark in addition to 
the valence one. 

2.3 Structure function 

Now we have all necessary to calculate the pion structure function ^(x). This function is given by the con- 
volution of the pion light-cone soft wave function and the corresponding constituent quark structure function 

KM 



F 2 (x,Q 2 )= Yl C dkldx'^(x\k ± ) 2 fi(^,Q' 



(33) 



j=Q,Q 

The leading order expressions for these structure functions in terms of parton densities read 



F 2 (x)= &Qi(x), ti(x)= ]T e 2 xql{x) (34) 

flavours flavours 

where Qi and q\ are quark parton densities of a given flavour i inside a pion and constituent quark j, respectively. 

At large x the behavior of the structure functions is determined by the valence quark (antiquark) distribution 
v(x) inside the constituent quark (antiquark). The latter is taken in the form 

v(x)^x- a (l-x)- b (l-N{l-x) c ) (35) 

The parameters in this formula are chosen such that v(x) should satisfy the following sum rules at Q 2 , = 10 GeV 2 

l l 
J v{x,Q^)dx = l J v(x, Q^)xdx = 0.35 (36) 



The last number is the fraction of constituent quark momentum carried by the valence quark. Due to the sum 
rule for momentum, this value is equal to the fraction of hadron momentum carried by valence quarks. For a 
proton this fraction is known to be w 0.35 at 10 GeV 2 . 
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Table 1: Table of the parameters in the pion form factor calculations. We assumed tog(k 2 ) = m G (l — k 2 /0.5) 
as k 2 < 0.5 GeV 2 and ma — as n 2 > 0.5 GeV 2 , where k is the gluon momentum. 





Set 1 


Set 2 


Set 3 


Set4 


constituent quark mass m, GeV 


0.35 


0.35 


0.35 


0.35 


m G , GeV 


0.7 


0.7 








a s in (28) 


a s (Q74) 


a s (n 2 ) 




a s (K 2 ) 


9a(0) 


0.86 


0.86 


0.86 


0.86 



3 The results of fitting the data on form factor and deep inelastic 
parton distribution 

Fig.| shows the results of fitting the data on F„(Q 2 ) |(J via the formulas @ and <^) for F ss and for F sh . 
Fig. 3 a demonstrates the region of low Q 2 where F ss dominates; Fig.||b emphasizes the region of moderately 
large Q 2 , Q 2 = 5 -j- 20 GeV 2 . Fig.^ plots the typical contribution of the soft-hard term to the pion form factor. 
Several variants of the calculation relate to different sets of the parameters tuq , a s , and m presented in the 
Table |l|. Let us point out that g A calculated with the determined wave function appreciably depends on the 
constituent quark mass: for a heavy constituent quark g A is rather small and equal to 0.86, whereas for a light 
constituent quark g A is close to unity. 

In all the fits the soft pion wave function ^{x^k 2 ^) was the basic variational quantity. Fig. || presents the 
typical reconstructed soft wave function. A specific feature of the reconstructed wave function is a quasi-zone 
character of gg-excitations in the pion. We tried to find a parametrization without a dip in the region 2-^4 GeV 2 
but failed; all the used fitting procedures suggested qualitatively similar double-humped wave functions. So 
one can think this dip to reflect some essential feature of the qq dynamics inside the pion. However, we have 
no definite ideas on the origin of such a specific behavior. 

Fig.^j displays the distribution amplitude (f>(x) calculated with the determined G v (s) . It turned out to be 
very close to the asymptotic function (Q). 

Fig.^ presents the valence quark distribution inside the pion. The parameters of the distribution ( |35| ) b = 0.99 
and a and c from the range a = 0.5-^0.7 and c = 1-^2 were found to provide a reasonable description of the data 
Q. Mention that b is very close to unity. If the valence quark distribution is determined by the Reggeized- 
gluon exchange, then our result gives for its intercept the value close to unity in qualitative agreement with 



4 Conclusion 

We have analyzed the pion structure within a dispersion relation formulation of the light -cone technique when 
the pion in the soft region is treated as a two constituent quark bound state. At small values of the qq light-cone 
energy the pion is described by a model soft wave function, whereas at large values the one-gluon hard exchange 
is taken into account. This provides the form factor high-Q 2 asymptotic behavior in agreement with pQCD. 
The obtained pion form factor describes well the available experimental data at Q 2 10 GeV 2 . 
Our results are as follows: 

• We considered the soft-soft (F ss ) and soft-hard (2F sh ) contributions to the pion form factor within the 
light-cone quantum mechanics. The derived expressions involve the soft radial wave function of the pion 
which has been treated as a variational parameter of the approach. By fitting the data on the pion form 
factor at Q 2 < 3 GeV 2 , we determined this soft radial wave function. This allowed a calculation of 
the relative soft-hard contribution to the pion form factor in a broad range of momentum transfers. It 
turned out to be relatively small (less than 50% at Q 2 — 20 GeV 2 ) because we used a large value for the 
boundary of light-cone energy squared dividing the soft and the hard regions (s = 9 GeV 2 ) and hence 
we related a large portion of the pion form factor to the soft-soft contribution. However, smaller values 
of this boundary do not change qualitatively the results, except for quantitative increasing the soft-hard 
fraction. 

The calculated pion axial-vector decay constant agrees well with the experimental value. 

• The soft radial light-cone wave function as a function of the square of the light-cone qq energy s has been 
found to demonstrate a specific behavior: it is large at s < 2 GeV 2 , close to zero as s = 2 -f- 4 GeV 2 , and 
has a bump in the region s = 4 -j- 9 GeV 2 . Our attempts to find a wave function of a more regular shape 



8 




Figure 3: The pion form factor. Solid line - calculations with the Set 1 of the parameters; dashed-dotted line 
- Set 2, dashed line - Set 3, dotted line - Set 4. 
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Figure 4: The contribution of the soft-hard term to the pion form factor. The curve notation is the same as for 
the form factor. 



10 




Figure 5: The soft light cone wave function of the pion. 
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Figure 6: The distribution amplitude. Solid line - our calculation; dashed line - asymptotical curve (f>(x) = 
Qx(l — x). 
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Figure 7: Valence distribution in a pion xV(x, Q 2 = 10 GeV 2 ). The solid curve - a = 0.69, b = 0.99, c = 1; 
the dashed curve - a = 0.46, 6 = 0.99, c = 1.55. 
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failed as all the fits suggested a double-humped behavior. We have no definite ideas about the origin of 
such a quasi-zone character of the qq excitations. Nevertheless, such an unexpected wave function seems 
to reflect some unknown essential details of the quark dynamics in the pion. 

• The distribution amplitude 4>{x) calculated with the obtained soft wave function turns out to be very close 
to the asymptotic function 4>^ s (x) predicted by pQCD. 

• Describing the data on deep inelastic scattering off the pion allowed investing the parton structure of the 
constituent quark. The distribution of the valence quark-parton in the constituent quark is found to be 
in a qualitative agreement with the parametrization suggested by the Reggeized-gluon exchange. 

The authors are grateful to the International Science Foundation for financial support under grant R1000. 



5 Appendix A: Bound state description within dispersion relations 

To illustrate main points of the dispersion approach we consider the case of two spinless constituents with the 
masses m interacting via exchanges of a meson with the mass p. We start with the scattering amplitude 

A(s,t)=<k[,k' 2 \S\ki,k 2 >, s=(k 1 + k 2 ) 2 , t={k 1 -k[) 2 (37) 

The amplitude as a function of s has the threshold singularities in the complex s-plane connected with elastic 
rescatterings of the constituents and production of new mesons at 

s = Am 2 , (2m + /i) 2 , (2m + 2^) 2 . . . (38) 

We assume that an S- wave bound state with the mass M < 2m exists, then the partial amplitude A (s) has a 
pole at s = M 2 . The amplitude A(s,t) has also i-channel singularities at t = (np) 2 ; n = 1,2,3. . . connected 
with meson exchanges. If one needs to construct the amplitude in the low-energy region s > 4m 2 the dispersion 
N/D representation turns out to be convenient. Consider the S'-wave partial amplitude 

l 

A (s)= dzA(s,t(s,z)), (39) 



where t(z) = — 2(s/4 — m 2 )(l — z), z = cos9 in the c.m.s. The Aq(s) as a function of complex s has the 
right-hand singularities related to s-channel singularities of A(s,t). In addition, it has left-hand singularities 
located at s — Am 2 — (np) 2 ; n = 1,2,3.... They come from ^-channel singularities of A(s, t). The unitarity 
condition in the region s ~ Am 2 reads 



1 / 4m ^ 

lmA (s) = p(s) \A ( S )\ 2 , p(s) = —^1-— (40) 

with p(s) the two-particle phase space. The N/D method represents the partial amplitude as A (s) = 
N(s)/ D(s), where the function N has only left-hand singularities and D has only right-hand ones. The unitarity 
condition yields 

oo 

f ds p(~s)N{~s) 

D(s) = l- / — _ = l-B(s). (41) 

J IT S — S 

4m 2 

Assuming the function N to be positive we introduce G(s) = \J N(s). Then the partial amplitude takes the 
form 

Ms) = G(-s) [1 + B(s) + B 2 (s) + B 3 (s) + ...} G(s) = ■ (42) 

This expression can be interpreted as a series of loop diagrams of Figj^ with the basic loop diagram 

oo 

J 7T S — S 

4 m 2 
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Figure 8: One of the terms in the expansion of Aq(s) 




Figure 9: One of the terms in the series for T^. 

The bound state with the mass M relates to a pole both in the total and partial amplitudes at s = M 2 so 
B(M 2 ) = 1. Near the pole one has for the total amplitude 

A = < k[,k' 2 \P > - r0 1 < P\k u k 2 > +rcgular terms 
M — r 

= X* P (K,k' 2 ) M2 1 _ p2 Xp(k 1 ,k 2 ) + ... (44) 

where Xp(ki, k 2 ) is the amputated Bethe-Salpeter amplitude of the bound state. The dispersion amplitude near 
the pole reads 

A = N/D + regular terms related to other partial waves 

G 2 (M 2 ) G 2 {M 2 ) 
(M 2 - s)B'{M 2 ) M 2 - s ( 0) 

where G v is a vertex of the bound state transition to the constituents. The singular terms correspond to each 
other and hence 

Xp( fcl ,fc 2 )^G„(P 2 )^-^= (46) 

Underline that among right-hand singularities the constructed dispersion amplitude takes into account only the 
two-particle cut. 

Let us turn to the interaction of the two-constituent system with an external electromagnetic held. The 
amplitude of this process T M =< k[, k' 2 \J fJ- (q)\ki, k 2 > in the case of a bound state takes the form 

T, = <k[,k' 2 \P'> pl2 \ M2 <P'\Mq)\P> p2 ^ M2 <P\h,k 2 >+... 

= X*p(k[,k' 2 ) — 2 * w (P' + P),F(q 2 ) p2 _ X — XP ( kl ,k 2 ) + ... (47) 

where the bound state form factor is defined as 

<P'\J^q)\P>^(P' + P) fl F(q 2 ) (48) 

The dispersion amplitude T M with only two-particle singularities in the P 2 - and P /2 -channels taken into 
account is given 23 by the series of graphs in Figj9[ 



These graphs are obtained from the dispersion scattering amplitude series by inserting a photon line into 
constituent lines. The amplitude reads 

T,{P', P, q) = 2P t2 (q)T(s', s, q 2 ) + ^C, (49) 
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Figure 10: Pionic dispersion loop graph B W (P 2 ). 



P 2 = 8, P' 2 = S', q = P'-P, PM = ( P - S f ff)/i 

The dispersion method allows one to determine T(s, s', q 2 ), which is the part of the amplitude transverse with 
respect to q^. Summing up the series of dispersion graphs in Fig.^j gives 

T{S ' q } = — B(sj r(s ' S ' q ] 1^W)- (50) 

Here 



7r(s — S) 7r(s' — s) 

and A(s', s, q 2 ) is the double spectral density of the three-point Feynman graph with a pointlike vertex of the 
constituent interaction. 

The longitudinal part C is given by the Ward identity 



G(s) D ,„„ G(s') 

l-B(s) 



C= T ^(BW-B{s)) T ^ r (51) 



At s = s' = M 2 , the quantity develops both s and s' poles, so 

T,(P>, P, q) = ^t^{P' + P),F(q 2 )^^ + less singular terms (52) 

where 

2 f dsG v (s) dsG v (s) , 2 
F{q )= ] ^)^tf) A(S ' M } - (53) 

is the bound-state form factor (see (Q) and (]47|)). So, the quantity < P'| J M (<7)|P > corresponds to the three- 
point dispersion graph with the vertices G v . The following relation is valid A(s',s,0) = n5(s' — s)p(s). This 
is a consequence of the Ward identity which relates the three-point graph at zero momentum transfer to the 
loop graph. This relation yields the charge normalization P(0) = 1. The expression ( |53| ) gives the form factor 
in terms of the iV-function of the constituent scattering amplitude and double spectral density of the Feynman 
graph. In general, the following prescription works: to obtain the dispersion expression spectral density in 
channels corresponding to a bound state, one should calculate the related Feynman graph spectral density and 
multiply it by G v . 

Mention that only on-shell constituents contribute to all the quantities. If the constituent is a nonpoint 
particle, the expression (p3[) should be multiplied by form factor of an on-shell constituent. 



6 Appendix B: Pion within the dispersion approach 

We start with describing the quark-pion vertex. For on-shell quarks there is one independent structure 
iq a lbq a I ' with a a color index and N c = 3 the number of quark colors. We consider a ir + and omit 
the flavor which gives the unity factor. Let us introduce the momentum P^, P 2 = mj. We shall also use an 
off-shell pion momentum P, P 2 = s. The pionic dispersion loop graph Figjl^ reads 

B^P 2 )- J t G2 p 2 W ^ S *K) = 1 (54) 

J 7r(s - P 2 ) 

4m 2 
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Figure 11: The soft-soft contribution F£ s (q 2 ) to pion form factor. 

with p„(s) the spectral density of the Feynman loop graph 

P-n{s) = — ^— j j dkidk 2 5{k 2 — m 2 )5(k 2 — m 2 )S(P — fa — k 2 ) Sp ((fa + m)i-f 5 (m — k 2 )ij5^j 



s I 4m 2 „ 



= —\l 6{s-Am z ) 

87rV s k ' 

with m a constituent quark mass. 

The form factor of a pion is given by the following matrix element 

< Pi\J™(0)\P„ >= (Pi+P^ F^q 2 ) 

P% = P? = M 2 , (P^-P v ) 2 =q 2 . 
The double dispersion representation for the form factor corresponds to Fig.[ll| 

2 f dsG v (s) d S 'G v {s') nf.\-. G ( s ) 



1 

8^ 



Here f c is a constituent form factor, / c (0) = 1, and is defined as 

dfadk^dk^kl - m 2 )S(k' 2 - m 2 )8{kl - m 2 )8{P -fa- k 2 )S{P' - k[ - k 2 ) 

x Sp ((k[ + m)7 Al (fci + m)ij 5 (m - k 2 )ij 5 ^j = 2P^(q) A 7r (s', s, q 2 ) 

with 

The trace reads 

-Sp ^(fci + m)7 M (fci + m)7 5 (m - ^2)75) = ™ 2 (fci + + k lf i(k[^) + k'^fakz) - k2 fl (qfa 
Multiplying both sides of ( p8| ) by P M and using ( |59| ) one obtains 

with A(s', s, q 2 ) = (s' + s - q 2 ) 2 - 4s' s. At q 2 = one finds 



A 7r (s',s,g 2 = 0) = -\/ 1 S(s' - s) = np^(s) S(s' - s), 



and 

ds G 2 v {s) 
7r(s — m 2 ) 

As we have mentioned this is just the Ward identity consequence 



F ff (0)=/ C (0) / v l 2 \ 2 PAs) = l 
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One can equivalently formulate the dispersion approach on the light-cone by introducing the light-cone 
variables 

k- = -j=(k - k z ); k+ = ^=(fc + *,); k 2 = 2k + k- - k]_; (63) 

into the integral representation for the form factor spectral density ( |58| ) . Performing fc_ integration and setting 
(/i = +) in both sides of ( |59| ) one finds 

a / / i\ 1 f dxd 2 k± ( m 2 + k 2 \ ( . m 2 + (k± - xq±) 2 \ .. , , 2 . 

A w (s', s, q 2 ) = — — -6 s - — i Us' \- {(s' + s)(l - a) + g 2 a;) (64) 



167r j cc(l — x) 2 \ x(l — x) J \ £c(l — x) 

Here we denoted a; — k 2 +/P+ and fcj^ = k 2 ±. 

Substituting (Q) into ( |57| ) and performing s and s' integrations, one derives 

K(<li) = - I dxd 2 k ± tp(x,k ± )ijj(x,k ± - xq ± )fi(x,k ± ,q J _)f c (q 2 _), (65) 



7T . 

where the radial light-cone wave function of a pion is introduced 

?/>(£, fc±) = — = , s = — -r- ^ (66) 

(3= - ^ 2 + fc i-^ fc = 0) = l 

The quantity /3 accounts for the contribution of spins. It is different from unity because both the spin-nonflip 
and spin-flip amplitudes of the interacting quark contribute. The eq.(|62|) is the normalization condition 

dx dk\ \ip(x,k ± )\ 2 = 1. (67) 

Let us now consider the pion axial- vector decay constant f^. It is related to the 7r — > [iv decay as 

< 7 r+(F x )|A; a (0)|0>-i(P.) Al ^, U = 130MeV. (68) 

To derive the expression for this matrix element we must first consider the quantity 

< UD\A™{0)\0 > (69) 

with U and D constituent quarks and then single out the pole corresponding to the pion. Mention that the 
axial current A M (0) = ^(0)7^75(^(0) is defined through current quarks. 
The bare matrix element reads (Q denotes a constituent quark) 

< Q(k)Q(P - fc)|^ a (0)|0 > bare = Q(P - k) [W59° A (P 2 ) + P^ 5 h° A (P 2 )] Q{k) (70) 

If current quarks were identical to constituent ones we would have had 

g° A (P 2 ) = l, h A (P 2 )^0. 

It is reasonable to assume that at least g A (0) and h A (0) are not far from these values. The bare matrix element 
enters into a single loop graph whose spectral density reads 



/Nr. 



dkidk 2 6(k 2 - m 2 )5{kl - m 2 )6(P - ki - k 2 ) (71) 



xSp (hnJ 5 g A (P 2 ) + P»lsh° A (P 2 )}(h + m) i 75 (m - k 2 ) 
So the loop graph is equal to 

B, = iP»[Am 9 \{P*) - ^h A (P 2 )]^N c J ^^r^Vl^n (72) 

After allowing for constituent quark rescatterings we come to the series of dispersion graphs of Fig.^j which 
gives 
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B^ 




B 



B,, 



-+... 



Figure 12: The series of dispersion graphs for < QQ\A^(0)\0 >. 



< Q(k)Q(P - fc)|^ a (0)|0 >= Q(P - k) [ 7m 753a(P 2 ) + P^h A {P 2 )} Q(k) 



with 



g A (P 2 ) = g° A (P 2 



ds G(s) 1 



MO- -T^)M 9 5(P 2 )-£^(P 2 )]/- (8 _ M2)16 , 

The form factor /i^ develops a pole at P 2 = m 2 as P^rn 2 ) = 1. Near P 2 = to 2 one has 



-y/l - 4m 2 /' 



< QQ|A M |0 >=< QQ|tt > 



— p2 <7T|^|0> 



Comparing the pole terms in (34) and (35) and using the relation 

< n\QQ >= ^fic v 



f N e . 



one finds 



< 7r|^ a (0)|0 >= iP^m[g A {mi) - -^h A (mi)W N t 



2m 



and hence 



with 



J ir(s — M z ) 1d7t 



K = 9\{rnl)-'^h A {ml) 
2m 



(73) 



(74) 



7T(S — M- 2 ) lD7T 



(76) 



(77) 



We can neglect the second term because the small value h A is further suppressed by m 2 . Finally, one finds 

«« ff °(m 2 )« g° A (0)« 0.75- 1 (78) 
In terms of the light-cone wave function ([76]) takes the form 



V2tt 



dec dk i 



^/to 2 + fc 2 



ip(x,k±) = U 



(79) 



7 Appendix C: Soft— hard form factor 



The soft-hard form factor F sh , given by the graph of Fig.|l3|, accounts for the assumption that at s" > sq the 
QQ interaction is described by the convolution of the one-gluon exchange kernel with the soft-region vertex 
G v (the right block in Fig.^3j) . To derive the spectral density of F sh , we start with the corresponding Feynman 
graph with pointlike vertices 



2P„(«) F*V) = ^%^9 2 



dk dk' 
(2tt)H (2tt)H 

1 



(80) 



m 2 — k 2 m 2 — (P — k) 2 m 2 — (P' — k) 2 m 2 — k' 2 m 2 — (P' — k') 2 \i 2 — [k' — k) 2 
xSp ((kx +q + m)7 Al (fci + m)ry 5 (m - k 2 )j a (m - fc^Hs^'i + rn)l c 
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Figure 13: The soft-hard contribution F£ to pion form factor. 

To allow only for two-particle intermediate states we consider the contribution of s, s", and s' cuts and neglect 
the contribution of three-particle intermediate states with cutting the gluon line: three-particle states are 
beyond the scope of our consideration. The resulting spectral density over both soft s and s' reads 

f Hi" 

A sh (s', s, q 2 ) = 9{s < soW < s ) g 2 C F \ ¥ 6(s" > s Q ) Tr(s, s', s", q 2 ) (81) 



x 



dfadk 2 dk 1 . 2 2,.-, 1.2 2,w;//2 2 



8tt 2 



5{k{ - m 2 )S(k 2 - m 2 )5(k'{ 2 - m 2 )6(P -fa- k 2 )S(P" - k" - k 2 )S(k'{ -fa-q) 

f dk' dk' 1 

x / ^f^8{kf - m 2 )S(k' 2 - m 2 )6(P' -fa- k' 2 )- 



8^ 2 V1 JVZ ,y 1 2J m 2 G -{k 2 -k' 2 Y 

with P 2 = s, P" 2 = s", P' 2 = s'. 

The quantity Tr(s, s', s" ', q 2 ) is the dispersion expression for the fermion-loop trace with all fermions taken 
on mass shell (for details see the Appendix D) 

Sp ({fa +q + m)7 M (£i + m)i l5 {m - k 2 ) la {m - k' 2 )ij 5 (k[ + m) 7 Q ) = 2P^{q)Tr (82) 

Tr = 2s' (a(s", s, q 2 )(s" + a - q 2 ) + q 2 ) - 4m 2 («(*", s, q 2 )(s' + s - q 2 ) + q 2 )) 

where 

-q 2 (s" + s-q 2 ) 



ot(s",s,q ) 



(s" -s) 2 -2q 2 (s" + s) + q 4 

Use again the light-cone variables JlS)). Performing fc_ integration and denoting k = &2, k' = k' 2 , x = k 2+ /P + , 
x' = k' 2+ /P +1 k± = k 2 ±, k' ± = k' 2± we come to the final expression 



F-M) = 4*C F { ^ G v (s)0(s<s ) ;K> ? ) (83) 

' x(l-x) ~ m * x(l-x) 



where 



and 



dx'd 2 k' ± G v (s')8(s' < s ) a s (-(k-k') 2 ) 
16ttV(1 - x') m'+(fci-x' gx ) a _ 2 { ' ' ' q ' m 2 - (k' - k) 2 

X'il—x') 

,, 2 _ (k' ± x - k±x') 2 + m 2 (x - x') 2 
— {k — k) — 



m 2 + (k± — xq±) 2 , m 2 + (k' ± — x'q±) 



s" = r-— , s = 



x(l — x) ' x(l — x) x'(l — x') 

In ( |83] ) the renormalization g 2 — > 47ra s (— (k — k') 2 ) is taken into account. 

The distribution amplitude which describes the large-Q 2 behavior of the form factor (see |3l]) is expressed 
through G v as 
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Figure 14: Two loops in F sh : momentum notation. 



8 Appendix D: Trace calculation in the soft— hard contribution 



The soft-hard contribution is described by two-loop graph, and the dispersion technique prescribes that the 
total momenta squared of the QQ pair should be taken different for each loop, and then the integration over 
these values should be performed. So, we must first make the Fierz rearrangements to obtain trace calculations 
related to different loops. Namely, we group the expressions as follows 

2P fl (q)Tr = Sp ( (k'{ + m)7 (tl (fci + m)ij 5 (m - k 2 )j a (m - ^3)175(^1 + m)j a 



E 



Sp (h" + m)7 (tt (fei +m)«75(m- fe) la {m - k' 2 )iih(k\ + rn) 
G, x Sp 
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(k'{ + m)^^(ki + m)i-j^(m ~ k 2 )Oi Sp (m — /c^TsC^i + m)Oi 

i=S,V,T,A,P 

with Cs = 1, Cy = — 5, Ca = h, Cp = — 1. The second trace is nonzero only for A and P and we find 
2P /1 (q)Tr = ^Sp (hi + q + 01)7^ (fa + 771)175 (m - k 2 )l5la Sp (m - k' 2 )i^ 5 (k[ + m)^ a j 5 



-Sp 



(hi +g + m)7 ft (fci +m)i7 5 (m- fc 2 )75 (m - fe^Hs^i + m )7s = 2P /i (g)(Tr A + Tr P ) 



Each of the expressions Tri (i — A, P) is represented as a product of two factors, related to two different loops 
(see Fig.HJ) We must use the following relations 



for the left loop and 



hi + k 2 =P,P 2 = s, k'[ =h+q, k'[ + k 2 = P", P 
k{ + k 2 =P' \P' 2 = s' 



II D //2 _ // 



for the right one and set all the fermions on mass shell. This procedure yields 

Tr A = -4m 2 [a(s", s, q 2 )(s' + s - q 2 ) + q 2 } 
Tr P = 2s'[a{s", S, q 2 ){s" + s - q 2 ) + q 2 } 

And the final result reads 

Tr = [2s' (a(s", s, q 2 )(s" + s - q 2 ) + q 2 ) - Am 2 (a(s" , s, q 2 ){s' + s - q 2 ) + q 2 ))] 
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